General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA Technical Memorandum 78708 


HIGH FREQUENCY SOUND ATTENUATION IN 
SHORT FLOW DUCTS 

(NASA-TH-78708) HIGH FBECHENCY SOHND N78-23876 

ATTENDATICU IN SHCPT FLOW DOCTS (NASA) 25 p 
HC A 02/NF A 0 1 CSCI 20A 

Hnclas 

G3/71 16770 


J. W, Posey 


May 1978 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton'Virgima 23665 




SUMMARY 


A geometrical acoustics approach is proposed as a practical design 
tool for absorbent liners in such short flow ducts as may be found in 
turbofan engine nacelles. As an example, a detailed methodology is pre- 
sented for three different types of sources in a parallel plate duct con- 
taining uniform ambient flow. A plane wave whose wavefronts are not normal 
to the duct walls, an arbitrarily located point source, and a spatially 
harmonic line source are each considered. Optimal wall admittance distri- 
butions are found, and it is shown how to estimate the insertion loss for 
any admittance distribution. It is suggested that the extension of the 
methodology to realistic source distributions in variable area cylindrical 
or annular ducts containing arbitrary flow is conceptually straightforward 
and computationally practical on a vector-hardware digital computer. 

INTRODUCTION 

The study of sound transmission and attenuation in soft walled flow 
ducts has received considerable attention since the advent of jet-powered 
commercial transport aircraft in the late 1950's (see refs. 1 and 2 for 
reviews of this work). The foundation for much of this work was laid earlier 
by Cremer (ref. 3), who investigated the infinite uniform duct problem. 

Cremer concluded that regardless of the original pressure distribution of an 
acoustical disturbance, there will be a finite distance down the duct past 
which the lowest order, or least attenuated, mode will dominate the sound 


field. He suggested that overall sound attenuation would be maximized by 
max sizing the attenuation of the least attenuated mode. An excellent 
exposition and extension of Cromer's work is given by Tester in reference 4. 

Tyler and Sofrin (ref. 5) were the first to emphasize the importance 
of the higher order duct modes in the turbomachinery noise problem, and 
Zorumski (ref. 6) presented a systematic analytical approach to account for 
modal propagation (transmission and reflection) in finite ducts having 
multi-segment liners. Rice (ref. 7) has developed a liner design technique 
based on an observed collapse of infinite-duct modal characteristics with 
cutoff ratio, and Nayfeh, et al., (ref. 8) have derived a modal propagation 
algorithm for variable area ducts. 

As modal theory has been applied to more complex duct and flow con- 
figurations, physical int ?rpretations have become imprecise and computational 
problems (eigenvalue calculations most notably) have become more formidable. 
In attempts to circumvent some of the problems inherent in modal calcula- 
tions, a number of numerical field solution algorithms have been developed 
(refs. 9, 10, and 11). Such numerical procedures are powerful tools for 
attacking problems involving arbitrary duct shapes, arbitrary flows, compli- 
cated Sviu.ce distributions and variable boundary conditions. However, 
numerics go directly from governing equations to resultant pressure distri- 
butions, providing no direct physical insight into mechanisms of propagation, 
reflection and attenuation. 

Geometrical acoustics avoids many of the modal calculation proglems, 
and it also provides the physical insight lacking in purely numerical 
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approaches. Ray theory has been employed by Wright (ref. 12) to compare 
noise radiation from free field and ducted rotors, by Tester (ref. 13) to 
develop some insiqht into why modal attenuation is maximized in an infinite 
duct by use of the Cremer optimal impedance, and bv tester and Posev 
(ref. 13) to verify the optimal wall impedance and maxima attenuation 
calculated from modal theory for a point source in a finite cylindrical 
duct. These studies indicate that geometrical theories can be used with 
confidence in short ducts and when wavelengths are much less than the duct 
height or diameter, as is often the case in turbofan inlets and short 
exhaust ducts. As an example application, the present naper utilizes geo- 
metrical acoustics to develop design procedures for al""rptive liners in 
short parallel plate flow ducts for three types of high frequency acoustical 
sources: a plane wave with arbitrary angle of incidence on the walls, a 

point source and a line source. For each source type it is shown how to 
determine the wall admittance distribution which maximizes sound absorption 
and how to estimate the insertion loss for an arbitrary admittance distri- 
bution. 

PARALLEL PLATE DUCT 

The coordinate system for the parallel plate geometry is shown in 
figure 1. The two semi-infinite plates occupy the planes y = * h for 
z < d. There is a uniform flow in the z direction of Mach number M. 

Plane wave propagation . - First, consider the idealized case where the 
sound field at z = 0 can be represented as a plane wave with wavefronts 
normal to the y-z plane and making an angle 0 with the z axis, lhen the 


unit vector n normal to the wavefronts makes an angle P with the vertical. 
In the absence of flow, Huygen's Principle states that the wavefronts advance 
at the speed of sound c. In a uniform flow, the wavefronts still advance 
with the speed of sound in a reference frame moving with the medium, but the 
field is also convected with the medium at a velocity cMe, where e 2 is 
a unit vector in the z direction. Thus, the rays are no longer normal to 
the wavefronts, but rather propagate in the direction of the vector n ♦ Me z 
as shown in figure 2. The angle of incidence of the rays on the wall (the 
angle that the ray makes with vertical) is given by 3 R where 


tan P R 


M + sin P 
cos 9 


( 1 ) 


It is clear that only waves with P R in the range (0, it) carry energy 

toward the duct termination. Due to symmetry about 0 R = n/2 , attention is 

further limited in this section to P R in (0. ir/2) with no loss i nenerality. 

Therefore, tan P R > 0 or sin 0 > -M. If M < 0, then flow is coming into 

the duct and convection will prevent the escape of any acoustical energy 

unless the z component of n exceeds |M|. 

Tester (ref. 13) has considered this problem in some detail where the 

incident plane wave pressure field is given by p^exp(ik y v + ik z z - iwt) 

and the reflected wave by P r .exp(-ik y y + ik^z - iwt). Requiring continuity 

of particle displacement between thr wall and the fluid and defining the 

specific normal admittance in the usual manner, P = p c( ) , 

0 \ p 'wall 

Tester obtains a pressure reflection coefficient C r r p r /p- given by 
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( 2 ) 


2 

r = c os 8 - B(1 - M sin 0) 

p O 

cos 0 + 0(1 - M sin 0) 

Notice that complete absorption occurs when he admittance is 


6 op t = 0 - M sin 0)“ 2 cos 0 


(3) 


That is, the optimal duct liner has zero reactance and a specific resistance 
R opt equal to B'J t . 


R ODt = ^ M S ^ n 0 ) 2 /COS 0 


(4) 


Optimal resistance as a function of wavefront orientation is plotted 
in figure 3 for both inflow and exhaust. The Mach number magnitude is 
arbitrarily set to 0.35. These curves suggest that designing a liner for 
the attenuation of random incidence noise may be more easily accomplished in 
an exhaust duct than in an inlet duct due to the relative flatness of the 
curve for the exhaust case. That is, a normal specific resistance of unity 
would be near optimal for a wide range of wavefront orientations (-20° < 0 < 
70° in the present case) for propagation with the flow, while no resistance 
value would yield similar performance for propagation against the flow. 

For a specific design problem where a range of 0 or 0p is specified 
along with a range of Mach numbers, equations 1 and 4 can be used to con- 
struct a family of curves similar to those in figure 3 from which a target 
value for the resistance can be chosen. An estimate of the transmission 
loss for any given M, 0, duct height 2h and distance to duct termination d 
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is readily obtained by adding the fraction of the incident power which is 
radiated directly to the fraction which is radiated after a single reflec- 
tion, as illustrated in figure 4. Let be the fraction of the rays 

which radiate without reflection and the fraction which are reflected. 

r 

2 

The power reflection coefficient is |C r | , so that the insertion loss is 
approximately 

IL 5 -10 log(Fp + F r |C r | 2 ) (5) 

Multiple reflections can be similarly taken into account. 

A better estimate of the transmission loss could be found by employing 
ray scattering theory as presented by Felsen and Yee (ref. 15). 

Point source .- When a point source is radiating into a stationary medium, 
the resulting wavefronts are concentric spheres, as shown in figure 5a. If 
the source remains stationary, but the medium has a uniform flow Mach number 
M, then the wavefronts are still spherical, but they are no longer concentric, 
since they are convected downstream as they expand (see figure 5b). However, 
the ray paths are straight in either case (ref. 16). Thus, for a point source 
at (0, h s , 0) in the parallel plate duct of figure 1, the angle of incidence 
for direct rays on the duct wall at (x, ±h, z) is given by 

6 r * arctan (x !hVh )a 

s 

2 2 2 1 

The wavefronts may be considered locally plane when k x + z + (h ± h $ ) 2 

» 1, and an analysis along the lines of that used by Tester (ref. 13) in 
deriving equation 2 yields the following expression for the complex pressure 
reflection factor. 
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( 7 ) 


2 

cos 0^ - 6(1 - M cos 0 z ) 

C r ” cos 0^ +6(1 - M cos 0 Z ) 2 

Here, 0^ is the angle which the wavefront normal makes with the e y 
vector and 0 z is the angle it makes with e^. Equation 7 could be adjusted 

to account for the curvature of the wavefronts (see ref. 17), but this is 

beyond the scope of this paper. 

In order for equation 7 to be utilized in determining an optimal 
impedance distribution 8 0 p t (x, ±h, z), the wavefront orientation angles 
0^ and 0 z must be expressed as functions of position on the walls. It is 
convenient to define <t> as the angle which the ray makes with e z< That 
is, 

[(±h - h ) 2 + x 2 ] h 

$ = arctan (8) 


Since the ray propagates in a straight line with the velocity cn + cMe z> 
it can be shown that 


tan <t> = 


sin 0 z 
M + cos 0 


(9) 


Equations 8 and S combine to give 


sin 0 z 
M + cos 0 Z 



( 10 ) 


Since the flow is parallel to e , the center of curvature of any wavefront 
will always lie on the line y = h s » x = 0. Thus, once 0 z is determined 
from equation 10, 0^ can be found from 
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COS e y 


± sin 0 z 

"l + x 2 /(±h - h s ) Z ]* s 


( 11 ) 


Now that the wavefront orientation angles 6 y and can be determined 
for any wall position by using equations 10 and 11, equation 7 can be used 
to determine the reflection coefficient at any wall position for an arbitrary 
admittance B at that point. As with a plane wave source, can be 
nulled. In the present case, this is accomplished by 


cos 0 (x, ±h, z) 

B 0 p t (*. ±h, 2. M) =[,— M cos e'(x, *h, £T)2 ' 


( 12 ) 


Again, B Qpt is real, implying that the optimal reactance is zero and that 
the optimal resistance R Qpt is given by 8~ pt , which is now position 
dependent. 

If it is not possible tc va*"y the wall admittance continuously, as 
dictated by equation 12, then tne designer could divide the liner into uniform 
segments such that B Qpt had little variation over any one of tne segments. 

In determining the average value of B Qpt over a segment, the appropriate 
weighting function is the density of direct rays incident on the wall as a 
function of position (proportional to energy flux per unit frea, not to be 
confused with wave intensity). The ray density over one of the spherical 
wavefronts is inversely proportional to the square of the radius of curva- 
ture of that wavefront, but the surface of maximum density is locally normal 
to the direction of ray propagation. Let e r (x, y ; z) be a unit vector in 
the direction of the ray path, then the local density of rays th'ough the 
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surface normal to e^ is proportional to 
the density of rays incident on the wall is 




-1 


It follows that 


D(x, ±h, z) = A 


cos P R 



( 13 ) 


where A is the constant of proportionality and r c is the radius of curva- 
ture of the wavefront, given by 


r 

c 


[x 2 + (y - hj) 2 ]* 1 

sin 0^ 


( 14 ) 


Once a distribution of 3 has been selected by the designer, the result- 
ing insertion loss can be estimated by approximating P^, the power radiated 
from the lined duct, as the sum of the direct radiation and that carried by 
rays which escape the duct after one and only one reflection. Such singly 
reflected rays strike \.be walls at z > z t where 

, (h + h )d 

z ' = 2irnr ~ • y * ih < 15 > 


Thus , 


P L 5 


D(x, h, z) | C r (x , h, z) pdxdz + 


Z + lx, 


r 00 


D(x, -h, z) |C r (x, -h, z) | dxdz 


z" 


h r 00 


-h -« 


Ad f.2 , 2 , / 

d + X My 

r c (n • e r ) 


- h 


) 2 ]" li dxdy (16) 


9 


where C r is the pressure reflection coefficient given by equation 7. The 
inclusion of multiply-reflected rays is straightforward, but beyond the scope 
of this paper. The radiated power from the corresponding hardwalled duct is 


. ro-. .. .. + D(x, -h, z)j dxdz + 


J ) 

0 -OO 


and the insertion loss is 


h r® 


Ad 


-h -°° 


r‘ c Cn ■ e r ) 


[d 2 + x 2 + (y - h s ) 2 J“' s dxdy 


(17) 


IL = -10 log 



(18) 


Line source .- A line source is taken to lie parallel to e at z = 0, 

i (k x - ut) 

y = h.. The source flux varies as e s . Hence, in the absence of 

flow, the wavefronts are right circular cones and have a semi -apex angle 
given by arccos (k s /k), where k = u/c, and c is the ambient speed of sound. 
The ray density at any point on the wavefront is inversely proportional to 
the radius of the cross section of the cone, and the rays are normal to the 
wavefronts. Again imposing the uniform flow velocity cMe z , the wavefronts 
remain right circular cones, but the axes are no longer coincident with the 
source line. The intersection of a typical wavefront cone with the plane 
y = h $ is shown in figure 6. Huygen's Principle implies that the position of 
such a wavefront at any given time may be determined as the envelope of the 
spherical wavefronts (representing equal phase) from all the points on the 
source line which the wavefront has passed. Since these spheres are convected 
at the velocity cMe z as their radii increase at the rate c, and the apex of 
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the conical wavefront advances at a speed of co/k $ along the line source, 
it follows that the axis of the conical wavefront lies in the plane y = h t 
and makes an angle t with e x . 

tan t = Mk $ /k (19) 
Also, the semi -apex angle of the cone is y, where 


ta " Y ■ f 



(1 - M?) r 


•s 


( 20 ) 


Tnerefore, for k $ / 0 the instantaneous locations of the wavefronts are 
given mathematically by the family of surfaces satisfying the following 
equation. 


o 2 

(y - h )“ + [(a - x) sin t + z cos x] 

2 2 

= [(a - x) cos r - z sin xl tan Y, x <_ a (21 ) 

Here, (a, h s , 0) is the location of the apex of the cone. From equation 21 
it is straightforward to determine the relative transverse position of the 
apex (a - x) for wavefronts which are incident upon the walls y = ±h at 
any given axial station z. In particular, (a - x) is the positive root of 
the quadratic equation 

C 2 (a - x) 2 + C^a - x) + C Q = 0 (22) 
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where 


2 2 
C, = tan t - tan y 

= 2z tan t( 1 + tan^Y) 

(-h - h ) 2 2 22 

C = — ♦ z (1 - tan'i tan^Y) 

COS**T 

Let { be the length of the line between the point of incidence and 
the apex of the wavefront. Then from figure 6, it follows that 

t = [(a - x) 2 + (th - h $ ) 2 + z 2 ] S . (23) 

and cAt ■ t tan Y. Thus, the wavefront normal is in the direction of the 
vector f-om [x - i tan Y sin(t + Y), h $ , M itan Y] to (x, *h, z). 

n = ^ (24) 


where 


N - e i tan Y sin(T + Y) + e (th - h ) 
x y s 


+ e 7 |z - M ( tan Y] . 


and 


Since the components of any unit vector are its direction rosii.es, 

cos 6 = n • e 

y y 

cos 8 * n • e, (25)' 

z z 


12 



These cosine values can be substituted into equation 7, and nulling the 
numerator would yield an optimal axial variation of admittance for each wall 


e opt ( ±h » z » M > t ) ■ ( ±h - h s )n - H < z - M tan Y >i' 


(26) 


Notice that equation 26 differs from equation 12, the optimal admittance 
for a point source, in that there is no x dependence, but there is a depend- 
ence upon the ratio of the source wave number k s to the free space wave 
number k = Hence, it would be pasier to construct the optimal liner for 
the line source since it would need to vary in only one direction, but addi- 
tional information about the source (the values of w and k s ) would be 
required. 

The ray density on the wavefront is inversely proDortional to the radius 
of spreading, denoted by cAt in figure 6, and equal to 2. tan Y. As with 
the point source, the local ray direction must be determined in order to 
arrive at an expression for the density D L of direct rays incident on the 
wall. Specifically, e f is the unit vector in the direction of n + Me z , 
and 


D L (±h, z) 


B 

e.(±h, z)tan Y 



(27) 


Here, B is the constant of proportionality. Thus, if multiply-reflected 
rays are again neglected, the power radiated per unit length of the source 
from an off-optimal lined duct is 


P L = 


D L (h, z)|C r (h, z)| 2 dz + D L (-h, z) |C r (-h, z)| 2 dz 


' r h 




BU(y, d)tan y]" 1 


( 28 ) 
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and the power radiated per unit length from a hard walled duct is 


rd 


P H = 


[D(h, z) + D(-h, z)]dz + 


5U(y, d)tan Y]" 1 dy (29) 


The insertion loss is 


IL = -10 log -4 
P H 


(30) 


CONCLUDING REMARKS 

In aircraft turbofan engines, the noise that propagates through the 
inlet and fan exhaust ducts is due to several different mechanisms. Some of 
these mechanisms, such as the rotor interaction with small scale turbulence, 
or with a small scale inflow distortion, may be considered as point sources, 
while others, such as the rotor-locked pressure field, may be considered as 
a ring source or a series of ring sources. The liner design methodology 
developed above for point and line sources in a parallel plate duct can be 
straightforwardly extended in principle to point and ring sources in ducts of 
arbitrary shape and even to a random distribution of sources within such ducts. 
Nonlinear liner interactions, the effects of multiple reflections and scattering 
at duct discontinuities could also be included, and radiation patterns as 
well as insertion losses could be predicted. Although the implementation o* 
the methodology would involve extensive computer programing, the effects of 
basic acoustical phenomena such as refraction, diffraction, reflection and 
absorption could be easily followed through the calculations, a definite 
advantage over modal or purely numerical approaches. 
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Ray tracing solutions to such complicated problems as turbomachinery 
duct noise propagation have been avoided in the past possibly because of 
prohibitively large computation times being required; however, the recent 
advent of vector- hardware digital computers makes possible reductions in 
computation time for such problems of at least an order of magnitude. 
Therefore, it seems quite within reason that the principles of geometrical 
acoustics utilized in the example problem of this paper could be the basis of 
a practical design procedure which employs more realistic source, duct and 
flow models than are currently used in procedures based on modal theory. 
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Figure 2.- Vector relationship between the unit wavefront normal n and the 
associated ray direction for a plane wave in a homogeneous medium moving 
with uniform flow velocity cK£ . 




Figure 3.- Op 
of wavefront 








Figure 5.- Schematics of spherical wavefront spreading associated with a poin 
source: a) medium at rest, b) medium in uniform subsonic motion. 
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